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We discuss a non-perturbative renormalization of n-point Polyakov loop correlation 
functions by explicitly introducing a renormalization constant for the Polyakov 
loop operator on a lattice deduced from the short distance properties of 2-pint 
correlators. We calculate this constant for the 5(7(3) gauge theory. 



o 

CO 

f^i ■ 1. Need of renormalized Polyakov loop 

Universal properties of the confinement-deconfinement phase transition in 
QCD are often discussed in the quenched approximation (~ SU(N C = 3)) 
where the Polyakov loop L is treated as an order parameter. In a lattice 
regularized version, L(x) — n^ T =i ^( x o, x ), where U(xq,x) G SU(N c ) 
denotes the gauge link variable on a lattice of size N% x N T . The expectation 
value of the traced Polyakov loop is said to be related to the free energy 1 of 
^ ■ a single test-quark put into a gluonic heat bath (F q ) and that of the heat 

bath alone [Fq) via: 

W s <AT3 £^«> ^ y- = e-M-*>\ (frl ^ 1) (1) 

x 

Zg and Z denote the partition functions of the system with and without 
test-quark a at inverse temperature (5 — l/T. However, it is well-known 
that (L) contains (UV-) divergent self-energy contributions, i.e. in lattice 
perturbation theory, 

<L)~e- ClA ^ 2 +o(/) ; (2) 

where 1/aT = N T and c\ denotes the coefficient 2 of the expansion in the 
bare coupling g 2 . From this expression it follows, that taking the contin- 
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For a more formal definition of (L) and Z q see Ref. 5. 
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uum limit of (2) at fixed temperature leads to a vanishing Polyakov loop 
expectation value even in the deconfined phase: The divergence appearing 
on a lattice has to be removed by renormalization 3 . 

Why is a renormalization of (L) of interest? Firstly (L) acts as an 
order parameter for the confinement-deconfincmcnt phase transition which 
should have a physical meaning. Secondly a renormalized Polyakov loop 
can be used to construct effective actions for the hadron dynamics near the 
phase transition 4 . And thirdly Polyakov loop correlation functions refer 
to free energies containing the potential energy and entropy which are of 
interest in heavy quark physics. For example investigations of meson and 
hadron properties via the Polyakov loop and the Polyakov loop correlation 
functions on a lattice show the need for both, the renormalized Polyakov 
loop and the renormalized n-point Polyakov loop correlation functions. 

2. Towards renormalized n-point Polyakov loop correlations 

Recently we suggested 5 that a renormalization of the Polyakov loop can be 
obtained through the renormalization of the quark-anti-quark free energy 
calculated at short quark anti-quark distances. At short distances the finite 
temperature quark-anti-quark free energy is given by the zero temperature 
heavy-quark potential. From this property it follows that the divergent self 
energy contribution in the quark-anti-quark free energy at finite temper- 
ature can be removed through a matching of its short distance behaviour 
to that of the heavy-quark potential at zero temperature, which may be 
calculated perturbativcly 6 . Once having done so, also the large distance 
behaviour of the finite temperature free energy is fixed and can be used 
to extract the renormalized Polyakov loop as no additional divergences are 
introduced at finite temperature. Clearly to each renormalization prescrip- 
tion corresponds a renormalization constant (which in general depends on 
the number of colors and flavors) and it is worth noting that the Polyakov 
loop and the Polyakov loop correlation functions are composite operators 
and as such need their own renormalization constants and conditions. 

It turns out, however, that once having fixed the renormalization con- 
stant for the Polyakov loop operator, the renormalization scheme applies 
to any n-point Polyakov loop correlation function at finite temperature. 
This important feature is based on the (hierarchic) divergence structure of 
n-point Polyakov loop correlation functions which allows to remove the di- 
vergences through one single renormalization constant. It thus follows that 
it is sufficient to calculate only 2-point Polyakov loop correlation functions 
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at finite temperature in order to extract the renormalized Polyakov loop 
and the renormalization constant. 

In the renormalization scheme described above, however, the renormal- 
ized quantities are defined only up to an arbitrary constant (c) in the zero 
temperature potential. In terms of 2-point correlation functions we may fix 
the zero temperature potential relatively to the potential of the string pic- 
ture at large quark anti-quark separations, for instance we define c through 
V( C \(R ^> 1) ~ crR + c + 0(1/ R), where V( c )(i?) denotes the zero tempera- 
ture potential and a is the string tension. Usually one takes the standard 
Cornell potential (c = 0). 

In order to be definite, we define the renormalized Polyakov loop 
and its renormalization constant Z^(g 2 ) through 

Lf c) (x) = J] Z? c) (g*)U(x ,x) = {z^)) Nr L(x). (3) 

x = l 

A multiplicative introduction of the renormalization constant in the lattice 
operator respect the center symmetry of SU(N) and ensures that also the 
expectation value of the renormalized Polyakov loop can be used as an order 
parameter for the confincmcnt-dcconfinemcnt phase transition. Moreover, 
(3) implies 



(f,£«(x,) ... f,£» (x, w ) ... Tri?*(y w )) = (z£( S 2 )) "(-..) 



The behaviour of (4) at infinite quark separations can be used to define the 
renormalized Polyakov loop, for instance in terms of 2-point Polyakov loop 
correlation functions, 



poo (rp\ pR rn p\ 

ee ^lim ^ ! J =-ln|((zg )( . 9 2 )) L)| 2 = -ln|(Lf c) )| 2 .(5) 



We should note that (5) respects the usual colour structure of a quark anti- 
quark pair: It can form a colour singlet or colour octet state. However, the 
renormalized quark-anti-quark free energies coincide 5,7 in all colour chan- 
nels at infinite quark separation. Consequently, the renormalized Polyakov 
loop expectation value, which can be estimated from exp(— FF^/2T), does 
not reflect a colour structure and its magnitude is related to the free energy 
of a single heavy quark placed into a thermal gluonic heat bath. 
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Figure 1. The renormalized Polyakov loop expectation value as denned in (5) calculated 
on lattices of size 32 3 X 4, 8, 16. The matching to the zero temperature potential utilizes 
the standard Cornell form (c = 0). The inserted figure shows a comparison of (L R ) 
for the case c = (squares) and c = — T c (triangles) at TJ>6T C to (continuum) 1-loop 
perturbation theory. 



3. Discussion of lattice results in SU(N C = 3) 

Lattice results 5 of the renormalized Polyakov loop are shown in Fig. 1. The 
renormalized Polyakov loop expectation value acts as an order parameter 
for the confinement-deconfinement phase transition and is well behaved in 
the continuum limit. At high temperatures, (L R ) approaches values larger 
than unity, which is supposed to be its high temperature limit. However, 
the renormalized Polyakov loop depends (multiplicatively) on the fixing of 
the zero-T potential. In order to demonstrate this dependence and to study 
contact with perturbation theory 8 we have inserted data with different fix- 
ing constants, c = and c = — T c (see inserted figure). Although the renor- 
malization scheme depends on this fixing, the renormalized Polyakov loop 
expectation value approaches an 'universal' value at large temperatures: In- 
deed, for different c and c it follows liniT^oo ) = limT—>oo(LP c \). Since 
a shift of the zero temperature potential corresponds to a non-perturbative 
contribution in the renormalized Polyakov loop < L R ) we expect the renor- 
malized Polyakov loop to coincide with (continuum) perturbation theory 
at high temperatures where the shift does not influence the magnitude of 
{L R ). We therefore also expect that the renormalized Polyakov loop defined 
in our scheme approaches unity at infinite temperature. 
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Figure 2. The non-perturbativcly determined renormalization constant defined through 



The renormalization constant for the Polyakov loop can be estimated 
in terms of {{L R ) / (L)) 1 / Nt . We only note here, that the values of Z R 
(see Fig. 2) are independent of the lattice extent in time direction and 
compensate the divergence appearing in (2). For small couplings g 2 , f3 = 
6/g 2 6 [4.8,6.6], the data for Z H can be well described with a function 
f(g 2 ) = exp(Q 2 (8/3)g 2 + Q 4 .9 4 ) with Q 2 = 0.106 and Q 4 = -0.065. These 
values are of the same magnitude as they appear in lattice perturbation 
theory 2 . 

4. Conclusion 

We have shown that the non-perturbative renormalization scheme applies 
to n-point Polyakov loop correlation functions once having fixed the renor- 
malization constant for one of these correlation functions. This can, for 
instance, be done for the simplest Polyakov loop correlation function, the 
2-point function which describes the quark-anti-quark free energy. The 
renormalization of the free energies is equivalent to having fixed the cor- 
responding partition function from which other physical quantities like the 
potential energy 7 can now be deduced. 

References 

1. L. G. McLerran and B. Svetitsky, Phys. Rev. D24 (1981) 450. 

2. U. Heller and F. Karachi, Nucl. Phys. B251 (1985) 254. 

3. For a recent discussion and further references see also: P. de Forcrand and 
L.v. Smekal, Phys. Rev. D66-011504 (2002), hep-lat/0107018. 

4. A. Dumitru and R.D. Pisarski, Phys. Lett.B525 (2002) 95. 

5. O. Kaczmarek et al, Phys. Lett. B543 (2002) 41. 

6. S. Necco and R. Sommer, Nucl. Phys. B622 (2002) 328. 

7. F. Zantow et al, hep-lat/0301015; contribution to SEWM-2002 
(this conference). 

8. E. Gava and R. Jengo, Phys. Lett. B105 (1981) 285. 



